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Abstract 

We consider the general spin-1 SU(2) invariant Heisenberg model with a two- 
body interaction. A random loop model is introduced and relations to quantum 
spin systems is proved. Using this relation it is shown that for dimensions 3 and 
above Neel order occurs for a large range of values of the relative strength of the 
bilinear (-Ji) and biquadratic (-72) interaction terms. The proof uses the method 
of reflection positivity and infrared bounds. Links between spin correlations and 
loop correlations are proved. 


1 Introduction 

1.1 Historical Setting 

In this work properties of the spin-1 Heisenberg model are deduced using a random 
loop model hrst introduced in the work of Nachtergaele m- Random loop models 
have been around since the work of Toth ||2T1 and Aizenman and Nachtergaele HI. In 
El a lower bound was obtained on the pressure of the spin-j Heisenberg ferromag- 
net; this improved the bound of Conlon and Solovej 0. Sharp bounds have recently 
been found 0. The loop model presented in m applies to the spin-j Heisenberg an- 
tiferromagnet. Both spin models can be applied to higher spins; for a review of these 
models we refer, for example, to IfTTII . The work of Ueltschi Il2^ combined and ex¬ 
tended these loop models. It has recently seen attention for its usefulness in several 
aspects of quantum spin systems. In 12^ it is shown that there is long-range order in 
various spin systems, including nematic order in the spin-1 system. The work of Craw¬ 
ford, Ng and Starr ||2| on emptiness formation also makes use of the model, as does 
the work of Bjornberg and Ueltschi ||4] on decay of correlations in the presence of a 
transverse magnetic field. The loop model presented here comes from IfTSl . it is sim¬ 
ilar in flavour to the Aizenman-Nachtergaele-Toth-Ueltschi representation. See Refs. 
lU El EH [m [19] and references therein. 

Quantum spin systems are currently a very active area of research. The growth of 
popularity of probabilistic representations has allowed new methods to be applied to 
these systems with many interesting results. This work looks at the general SU(2) 
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invariant spin-1 Heisenberg model with a two-body interaction 

= - Z ■ ^y) + -^2 (s. ■ Sv)'). (1.1) 

{x,y}eE ^ 

Here we will have x e A cZ‘^ and £ the set of nearest neighbour edges. The operators 
S = {S^,S^,S^) are the spin-1 matrices, see section|^for details of the model. The 
work in ll^ shows that in the region 0 < J\ < \ J 2 the system exhibits nematic order 
in the thermodynamic limit if the temperature is low enough and the dimension is high 
enough. Nematic order was also shown independently using different methods in Il20l . 
It is also shown that if A is bipartite there will be Neel order for - 0 < J 2 low 
temperature. This corresponds to the occurrence of infinite loops in the related loop 
model. Alternatively in < 2 infinite loops should not occur, it is proved in ifT^ that 
this is the case for J 2 - 0, the extension to ^2 > 0 should be straightforward. The first 
proof of continuous symmetry breaking was shown by Frohlich, Simon and Spencer 
Ha for the classical Heisenberg ferromagnet (and hence antiferromagnet). This result 
was extended by Dyson, Lieb and Simon [Sl to the quantum antiferTomagnet. The 
result excluded the case d - 3 and S - j, it was extended to this case in the work 
of Kennedy, Lieb and Shastry ifTSll . These works all used the method of reflection 
positivity and infrared bounds. For information on reflection positivity see Refs. EE] 
[TOjini and references therein. The Heisenberg ferromagnet is not reflection positive 
and hence does not benefit from these methods. 


1.2 Main result 

In this article we use the method of reflection positivity and infrared bounds on a ran¬ 
dom loop model. Links between cotTelations in the spin model and probabilities of 
events in the loop model are also derived in section 5. We focus on the quadrant 
Ti < 0 < J 2 , see Fig. (in fact we need only work on the unit circle - I 

but the quadrant is more convenient pictorially). The following result concerning Neel 
order follows from Proposition |5.3| (a). Theorem |6.1| and the discussion that follows. 
For the precise statements see Section]^ 

Theorem. For A c a box of even side length and d > 3 there exists a — a{d) > 0 
such that for 7i < 0 < ^2 if—J\ IJ 2 > ce then there exists c — c(a, d) > 0 such that 

|AHoo A ^ 
xgA 

Furthermore a(d) 0 as d ^ 00 . 

Note that the “liminf” version of this theorem will be proved. The limits exist but this 
will not be proved here. It is shown in the discussion after Theorem |6.1| that this sum 
is positive if 

IdKd<(-4Ji)/{-Ji+4J2). (1.3) 

Id and Kd are integrals to be introduced in ( |6.3] l. Their values for various d are given in 
the table below. 
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Figure 1 : The phase diagram for the general SU(2) invaiiant spin-1 model. Regions that are shaded darker 
have rigorous proofs of the relevant phases. The line J\ < 0, /2 = 0 is the Heisenberg antiferromagnet where 
antiferromagnetic order has been proven f8|, Neel order extends into the dark yellow region. The dark blue 
region 0 < < ^J 2 has nematic order at low temperatures da, with Neel order on the line J 2 > 0, 7j = 0. 

The adjacent dark yellow region has been proved to exhibit nematic order in high enough dimension ED. 
Antiferromagnet order is expected here but is not yet proved. 


d 

Id 

Kd 

3 

0.349882 

1.15672 

4 

0.253950 

1.09441 

5 

0.206878 

1.06754 

6 

0.177716 

1.05274 


It can be shown iiiia that Id 0 and Kj ^ 1 as c/ —> oo and that both are decreasing 
in d. This means we can prove that the region where Neel order occurs will increase to 
the entire quadrant Ji < 0 < J 2 as d ^ oa i.e., the ratio a(d) is decreasing. In = 3 
there is Neel order at low temperature in the spin system for -J 1 IJ 2 < 0.46, this is a 
triangular region of angle 65° measured from the Ji axis. 

Reflection positivity for this quadrant is already known; for Ji < 0 = ^2 it was shown in 
in and for 7i = 0 < J 2 one can see, for example, M lemma 3.4 for an explicit proof. 
It was proved in jS) that Neel order occurs for Ji <0 - J 2 , it is clear the result extends 
to a neighbourhood of the axis J\ < 0 < J 2 with J 2 sufficiently small. However it is 
impossible to extend the result concerning Neel order any significant amount without 
some new results. This is where the loop model has been essential. Indeed in an 
infrared bound is obtained of the form 




1 


2{-Jx)e{k) 


(1.4) 


where (A, B)ouh is the Duhamel correlation function and e(k) - 2 SjLiCl-cos k,) fork € 
A*. Notice that this bound becomes weaker as l^il decreases (equivalently on the unit 
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Figure 2; Events of the process a) represents single bars, b) represents double 

bars and c) represents the uniform measure on vertical segments being either parallel 
or crossing. 


circle as I/ 1 I/I/ 2 I decreases). Transferring this bound to {SQSl)^(k) requires the Falk- 
Bruch inequality, which would involve dealing with the term ([5 [J 2 (S t ■ S 3 ,) ,5^]])/?. 


After some calculation one obtains correlations in Proposition 5.3 such as {SISySlSy)/^. 


Hence to work directly in the quantum system using the methods of IS] one must obtain 
good bounds on these correlations. Simple bounds such as taking the operator norm 
are not sufficient due to the weakening of ( |1.4[ l as I/ 1 I/I/ 2 I decreases. Without using the 
loop model it is not clear how to obtain such bounds currently. 


The random loop model is presented in sections]^ andThe spin-1 Heisenberg model 
is introduced in Sectionj^ In Sectionj^the connection between the loop model and the 
quantum system is proved. In particular it is shown how to write various correlation 
functions in terms of probabilities of events in the loop model; some of these correla¬ 
tions are also presented in Il22]| . In Sectionj^the main result concerning Neel order is 
presented and proved. 


2 The random loop model 

We now introduce the loop model presented in ifTSl . To begin we take a finite set of 
vertices. A, with a set of edges, & c {{ 2 i:,y)|x,y e A,x y). We associate to this lattice 
a new lattice. A, and edge set, £; 

A=Ax{0,1), (2.1) 

6 ^{{(x,i),(y,jWJ 6 (0, l),{x,y) 6 £). (2.2) 

There are two lattice sites in A for every site in A and four edge in £ for each edge in 
£. We will write xq, xi in place of (x, 0), {x, 1). 

For y 6 > 0 consider a process, py^ y^, consisting of a Poisson point process on £ x [0,y8] 
and a uniform measure on segments of A x [0,y6] between events of the Poisson point 
process. The Poisson point process has two events that we will refer to as ‘single 
bars’ and ‘double bars’. Note that this process is on the edge set £, the events define 
corresponding events on the edge set £. The single bars will occur at rate -27i and 
double bars at rate J 2 for Ti < 0 < /y- The rate for the single bars is written in this way 
to be consistent with the connection to the quantum spin system that will be introduced 
in Section The interval [0,/l] will be referred to as a time interval. The uniform 
measure is on two possibilities, “crossing” and “parallel”. How to build loops from 
these events is described in detail below, see Fig. [^for pictorial representations of the 
events. 
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Figure 3: A example realisation with loops coloured differently, here there are four 
sites in the underlying A and for this realisation |-C(a;)| = 3. 


We first define the single and double bars. Single bars occur at a point (x,y,t) for 
{.r,y) € &. We define the corresponding geometric event on £ as a bar joining xi 
and yo at time t. Double bars occur at a point (x, y, t) and the corresponding event on 
£ is a bar joining xi and yo and a bar joining ;ico and yi, both at time f. A loop of 
length [ is then a map j ; {0,pi'\per —> A x {0,fi\per such that y{s) + y(f) s + t, 
y is piecewise differentiable with derivate +1 where it exists. If s is a point of non¬ 
differentiability then 7 (i-i-)) e £. Loops with the same support and different 

parameterisations are identified. For a realisation w of pwe associate a set of loops 
as follows: Starting at a point (x,-, s) e Ax [0,^8] we move upwards (i.e. in direction of 
increasing s). If a bar is met at time t it is crossed and we then continue in the opposite 
direction from (yj, t), where yj is the other site associated to the bar. Each maximal 
vertical segment between bars (xo,xi) x [i, f] (i.e. bars involving the site x occur at 
times s and t and no bar involving x occurs for u such that s < u < t) is either parallel 
(nothing happens) or crossing (the sites xq and xi are exchanged). If time (3 is reached 
the periodic time conditions mean we continue in the same direction starting from time 
0. We denote by -C(w) the set of all loops associated to a realisation u). Loops are most 
easily understood pictorially, see Fig. Note that the loops could be defined via a 
Poisson point process on £ x [0,jS] where bars can occur between x,- and yj with each 
{i, i) being equally likely. However one would still need to introduce the crossing or 
parallel events so that it is still possible to have xq and xi in the same loop even when 
there is no bar occurring on any edge containing x. 

For this loop model we have partition function 

J (2.3) 

Here 0 > 0 is a parameter and pj^j^ is the probability measure corresponding to a 
Poisson point process of intensity —2J\ for single bars and J 2 for double bars. The 
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relevant probability measure is then 
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(2.4) 


We are interested in sets of realisations, oj, where certain points of A x [0,y6] are in 
the same loop. Probabilities of these events are connected to correlations in the spin-1 
quantum system presented in Section]^ they will be required in the proof of Neel order 
in Section]^ Particular events of interest will be denoted pictorially , see Fig. These 
events are defined and denoted as follows. 


a) The event that sites x, and y, are connected (in the same loop). Note that the proba¬ 

bility of X,' and yj being connected is independent of i and j. Denoted E{xi — yj\. 

b) The event that xq and xi are connected, yo and yi are connected but there is no 


connection from any x, to any yy Denoted E 


Hi- 


c) The event that xq and yo are connected, xi and y i are connected but xo and xi are not 


connected. Denoted E 


. We can also have xq and yi connected and xi and 


yo connected but xq and xi not connected and denote the event in the analogous 
way. These events both have the same probability. 


d) The event that all four sites xo,Xi,yo,yi are connected. Denoted E 


The definition of bars means that if a loop is followed starting from a point x, 6 A 
(by moving in either the up or down direction) then the direction it is travelling upon 
arriving at a point y^ e A in the same loop is determined only by the number of bars the 
loop has encountered between the sites. For example on a bipartite lattice defined by 
sublattices and A^ such that {x,y) e £ <=> x e J^A,y ^ Ag the direction that x, is 
left and yj is entered will be the same if x and y are in the same sublattice and different 
if they are in different sublattices. 


Sometimes the order in which sites are encountered along the loop will be important. 
In this case arrows will indicate the order that sites will be encountered in on following 


the loop (up to parameterisation). The events E 


f T 


, E 




are the 


events that all four sites are connected and are encountered along the loop in the order 


indicated by the arrows. For example the first event,^ 


t I 


, means that upon leaving 


site Xo if we encounter yo before encountering xi then we will then encounter yi and 
then xi before closing the loop. As this notation is potentially confusing (but also 
seemingly unavoidable) the reader will be told explicitly when the order is important. 
When wanting the probability of these events we will drop the E from the notation, as 
below. 


It is intuitively clear that P(xo—yo) decays exponentially fast with respect to ||x - y|| 

/'(K—'o\ /'^O ''o\ /-*0 ^o\ 

for /3 small. Hence P |^| and P 4 must also have exponential decay. P |hi-| 

—>1/ \J^l—>1/ \*l >1/ 

should depend weakly on ||x - y|| for small enough p. For ||x - y|| large enough the 
probability may approach P(xo—xi )^, it is not clear how to prove or disprove such a 
relation at this time. 
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) Xi 

Vj 

Xo Vo 

*>) HI- 

xo^yo xq — 7/0 

c) 4 d)\x\ 



xi yi 

xi-^yi xi - yi 


Figure 4; Pictures representing the set of realisations where the pictured connections 
are present. 

3 Space-time spin configurations 


In order to make the connection with spin systems we need the notion of a space-time 
spin configuration. The spin system we shall connect to is the spin-1 Heisenberg model, 
we shall make this connection via an intertwining that merges two spin-j models. For 
this reason we will take 6-2 from Section 2 (25 H-1 for 5 = 5 ). This is also the reason 
the lattice A has two sites for every site in AT It is also possible to represent the spin-5 
model for general 5 by merging 25 spin -5 models, this will mean A will have 25 sites 
for every site in A. See ifTSlI for more details. This generalisation together with some 
results analogous to the ones presented here should be straightforward once the spin -1 
model is understood. It is not immediately clear which results will still hold however, 
investigation is required. 


From now on we take the cubic lattice in Z'^ with side length L, denoted A/,, with peri¬ 
odic boundary conditions. The edge set, £/,, will consist of pairs of nearest neighbour 
lattice points. Precisely 



11-* “5ll = 1 or |xi - y,| = L - 1 for some / = 1, (3.2) 


Where ||x - y|| is the graph distance between x and y. A space-time spin configuration 
is a function 


O- : Ax [Q,P]per 



(3.3) 


cTxij is piecewise constant in t for any Xi. We further dehne Z to be the set of all such 
functions with a hnite number of discontinuities. For a realisation of the process u we 
consider cr that are constant on the vertical segments of each loop in £.{cJ) and that 
change value on crossing a bar. This restriction on conhgurations will allow to make 
the link with spin systems. We call such configurations compatible with at and denote 
by S*^'i(a;) the set of all compatible configurations. The following relation holds as 
we work on a bipartite lattice, meaning hxing a conhguration’s value at some (x,, t) 
determines the conhguration on the entire loop containing (x,-, f): 




(3.4) 


from which we can obtain 


f 


Pyi,y2(dw) ^ 1. 


(3.5) 


We further dehne the set Zd^(w) d compatible conhgurations along 

with conhgurations that hip spin at points (x,-, 0 ) or (yj, 0 ) (or both) but are otherwise 
compatible. 
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In the remainder of this section it will be convenient to ignore the condition that com¬ 
patible conhgurations flip value on crossing a bars as it adds unnecessary extra com¬ 
plication. Hence we further define to be configurations that are constant on 

loops (and hence do not flip value at bars). D denotes the set of con¬ 
figurations in along with configurations that flip spin at points (x,, 0) or (yj, 0) 

(or both) but are otherwise consistent with the definition of The reader should 

bear in mind in the sequel that the connection with the spin-1 quantum system requires 
spin flips at bars, with this in mind the modifications to the remainder of this section 
required to incorporate the spin flips are easy. 

As in we will later need a more general setting for the measure on space-time 
spin configurations. We consider a Poisson point process on £ x [0,/3] with events 
being specifications of the local spin configuration. We will consider discontinuities 
involving two pairs of sites (xo, xi,yo>yi)- The objects of the process will be a set of 
allowed configurations at these sites immediately before and after t. We can denote 
these events as 


o’xo.i-o'ijj- o-j,| 


(3.6) 


Implicit here is an ordering on A with x < y. An event A is a subset of {-1/2,1/2)^ 
and occurs with intensity l(A). More precisely we let t : !P({-l/2,1 jlf) R denote 
the intensities of the Poisson point process, denoted p,. Given realisation, of p, let 
2 (^) be the set of configurations compatible with ^ meaning that cr e 2(^) if 

^XQ,t+ f+ Cr^.Q f+ CTy^ 

6 A whenever f contains the event A at point (xq, xi , yo, yi , ?)> 

IT.vq,,- 0-j;j j- ,f- 

and cr^. , is otherwise constant in t. The measure is then given by with the counting 
measure on compatible configurations. We note that different intensities can give the 
same measure as in 12 ^ . for i and i' intensities it is shown in ^Q3\ that 

A(d^) r^ F(o-) = rA+,'(d^) ^ F(o-). (3.7) 

crests) 



We want to write the Poisson point process involving bars in terms of intensities of 
specifications of spins. We require that specifications corresponding to single and dou¬ 
ble sets of bars have intensity —2J\ and J 2 respectively. If we naively define T by 


( (a' a a b 
\iF~c c b 


-2Ju 


a a 


a a 


- Ji- 


(3.8) 


For any a, a', b,c, c' e (1/2, -1/2), where the first event corresponds to single bars and 
the second event to double bars. We see there is an overlap on the specification 

b a a b 
B~c c b 


so this assignment of intensities of specifications cannot be correct. Simply removing 
the overlapping case from one of the specifications will result in events not having the 
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required intensities. This suggests we should instead define i by 


a a 


a a 


= -27i, 


b a 

L ■ 


a a 


a b 
c b 


- Ji, 


(3.9) 


a b 


B~c c b 


= J2-2J1. 


For any a,a',b,c,c' e {1/2,-1/2). Now each specification is disjoint from the other 
two and single and double sets of bars have intensities -27i and J 2 respectively, as 
required. We also have t(A) = 0 for any other specification. Then 





(3.10) 


This representation can be used to show reflection positivity of the loop model, however 
in light of the famous paper of Dyson, Lieb and Simon IHl showing reflection positivity 
for Ji <Q - J 2 and the recent paper M showing reflection positivity for Ji - 0 < J 2 
we have reflection positivity for Ji < 0 < J 2 more directly, as we shall see. 


4 The general spin-1 SU(2) invariant Heisenberg model 

Let S e jN. A spin-5 model has local Hilbert spaces t/Tx = Observables 

are Hermitian matrices built from linear combinations of tensor products of operators 
on ®xeK^^x- Physically important observables are expressed in terms of spin matrices 
5',5^ and S^, operators on that are the generators of a (25H-l)-dimensional 

irreducible unitary representation of 5u(2). They satisfy the commutation relations 

= (4.1) 

y 

where a,f3,y e (1,2,3) and &cri 3 y is the Levi-Civita symbol. Denote S = (5',5^,5^), 
then S ■ S = 5(5 -t 1)1. The case 5 = | gives the Pauli spin matrices. For 5 = 1 there 
are several choices for spin matrices, we will use the following matrices: 


5 ' = ^ 

'0 

1 

o' 

52 = — 

'0 -i 

O' 


4 

0 

o' 

1 

0 

1 

i 0 

—i 

, 52 = 

0 

0 

0 

V2 

lo 

1 

Oj 

V2 


oJ 


lo 

0 -ij 


Consider a pair (A, £) of a lattice, A c and a set of edges, £, between points in A. 
We will take A to be a box in Z'^, hence A is bipartite. We denote by A^ and Ab the 
two disjoint lattices such that A^ U A^ = A and every e e S contains precisely one site 
from A^ and one site from Ab- 

Then we take the operator 5',. for i - 1,2,3 to be shorthand for the operator 5 (,®Wa\p). 
Recall the definition of A and £> above, we shall use these below. 

The most general SU(2) invariant Hamiltonian with two-body interactions for spin-1 is 

= - Z + -^2 (S. ■ S,)'). (4.3) 

{x,y}e6 ^ 
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We will soon drop the parameters J\,J 2 from for readability. In this article we 

will be concerned with the region where 7i < 0 < J 2 - Associated to this Hamiltonian 
we have the following partition function and Gibbs states for (3 > 0: 


A,/? 


-T re 




<•) 


J\ J2 
Afi 


1 


yJ\j2 

Afi 


Tr ■ e 


-HH] 


,J\ J2 


(4.4) 

(4.5) 


Again we shall drop the parameters Ji , J 2 from the notation. 

The following new definitions come from Nachtergaele US). We introduce an isometry 
y : ^ with the property VD^(g) - (D2(g))®^y for g € S U(2) and the 

spin-5 representation of SU(2). Here the representation D' is given by the matrices 
and is given by the Pauli matrices, it is clear such an isometry exists as we 
can define it for spin matrices and then extend by linearity (recall that the spin matrices 
generate the representation). From this we obtain the key relation 


VS‘ = (cr'® 1 H- 1 ®cr')y 


(4.6) 


where cr' are the spin-f matrices (hence 2cr‘ are the Pauli matrices). Further we have 


y*y = 1 and yy* = p. 


(4.7) 


where P is the projection onto the spin triplet. Hence y5' acts on ® and so using 
the notation before y^5 ^ acts on We make the following definition 

R' y5'y*. (4.8) 


One can check that R' = (cr' ® 1 - 1 - 1 ® cr'). To make expressions more concise we 
will also denote Ay 'Sxex^x for X c A. For these new operators we have a new 
Hamiltonian (note we have now dropped the Ji and J 2 parameters) 

- Z ■ ^y) + -^2 ^ (4-9) 

and associated Gibbs states 

^TrPAe^^^'^\ (4.10) 

(4.11) 

A^ 

The connection with the previous Gibbs state can easily be made explicit, 

(A) A./, = <VAAy*)^9 . (4.12) 

We use Dirac notation in the following way: \a, b) denotes an element of the one site 
Hilbert space ® and \a, b) ® |c, d) for two sites etc. 

There are two operators of particular interest, both act on two sites. Firstly we define 
by its matrix elements 

{a',b'\®{c',d'\S^^^a,b)®\c,d) = {-If-'’'6a,a'5d4'5b,-c6b'.^c'. (4.13) 
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Geometrically this requires spin b and c and the spins b' and c' to be the negative of 
each other and also requires a - a' and d-d'. This corresponds to the the single bars 
in the loop picture. The second operator, , is also defined via its matrix elements 

{a',b'\®{c',d'0^'>'\a,b}®\c,d} = (4.14) 

The geometrical interpretation this time is that of the double bars. The actual operators 
needed are P and = PD^^^ P in order to account for bars occurring 

between any Xj and yj with each i and j from {0,1) being equally likely. Note here that 
from this definition we see that we require the spin value to change sign on crossing 
a bar as was mentioned in Section |3 There are also extra factors in and of 
e'™ for the bottom half of a bar (denoted n) and e for the top half of a bar (denoted 
U) where a - ± j is the spin value on the site in associated to the bar. By direct 
computation of the matrix elements we can prove the relations 

= (4.15) 

-P,,y. (4.16) 

Using these relations we can rewrite the Hamiltonian in the region /i < 0 < 72 as 

- Z ■ (4.17) 

ix,y}e£ 

5 The random loop representation 

We can neglect the term (7i + J 2 )Px,y in the Hamiltonian ( |4.17| l and instead add (27i - 
72 ) 1 , this does not change the Gibbs states. Doing this allows to use a useful lemma 
from IT] 

expi- ^ {uAy,^y + vB^^y-u-v)\^ I p(dw)]~[ Qj,. (5.1) 

( {x,y}eE j (x,y)eoj 

Here p is the measure associated to a Poisson point process on fix [0,1] with two events 
occurring with intensities u and v respectively. The product is ordered according to the 
times at which the events occur. C is either A or B depending on which event occurs. 
This is actually a slight extension of the lemma presented in HI. From this we can 
obtain 


- Z + -^2^^ + 27i - 72 ) U r P(dcu) YY ^yj (5.2) 

( (A-,j>)e£ j ^ {Xi,yj)eb> 

here each A^' is one of or 2)^'. The process has intensity -27i for single bars and 
J 2 for double bars. Again the product is ordered by the time events occur. 

We now prove the connection between the loop model and the quantum system. This 
will enable us to understand certain important correlation functions. After this we 
should have the tools we need to calculate any two point correlation (at least ones in¬ 
volving only spin operators). The first thing to understand is the extra factor, which we 
shall denote by Zxijjicr, cu), the product of all factors e*'™ from operators and 2)^^ 
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corresponding to the bars in loop(s) containing x, and yj in a realisation a> of 
Again a e {1/2, -1/2) is the value that cr assigns to the portions of these loop(s) in the 
Ax sublattice (or if all of a loop is on the sublattice Ab a is given by the negative of 
the value assigned to the loop). The value of Zxi,yj(o', <A) is specified by the following 
lemma: 


Lemma 5.1. For A bipartite we have for all i, j 


Zxi,yj(a-, CJ) 


1 if a 6 

(_l)lh-rll if a- e \ and co e E[xi-yj\. 


(5.3) 


Before the proof we should note that the lemma says that the only dependence on cr is 
at Xi and yj at time zero. If the spin does not flip at both sites that we get total factor 
1, else it depends on which sublattices the sites are in. If the spin only flips at one site 
then there are no compatible configurations hence the value of the total extra factor is 
unimportant. 

Proof To begin note that we can take {i,f) - (0,0). The result for (/,/) (0,0) 

follows as the choice of i or j does not affect which sublattice the two sites are in. 
Suppose cr e E*''(fa;). Moving upwards from xq the first bar encountered is n, the 
bars encountered then alternate between U and n. Moving downwards from xo we first 
encounter a bar U then alternate between n and U. This means we can make a matching 
between bars of the form n and bars of the form U. Because there are no spin flips at 
time zero all the bars n have factors e'™ and all the bars U have factor where a is 
the spin value cr gives to xq at time zero. Hence we have full cancellation and are left 
with factor 1. If there were a spin flip then bars between xo at time 0- and yo at time 
0+ would have factors and for n and U respectively. 

If cr e E^*'j,^(w) \ 2^'^(w) and (-l)ll''~>’ll = 1 and to e ^[xo—yo], then xq and yo are 
in the same sublattice. We can thus deduce that the section of loop that moves up¬ 
wards/downwards from Xo crosses an even number of bars before reaching yo. This 
means that the loop containing xo and yo contains an even number of bars of each type 
(n or U). Hence we can make a matching of a bar U in one ‘half’ of the loop with a bar 
U in the other ‘half’ and the same with bars n, with some bars left over. The factors 
from bars in the matching will thus be 1 as the spin flip at xo at time 0 means one bar 
in each pair has factor and one bar has factor Here by ‘half’ of a loop we 

mean the section that connects xo at time 0+ with yo at time 0- or xo at time 0- with 
yo at time 0+. There are still possibly some bars left over as each half of the loop may 
have a different number of bars in it. A moments thought reveals that there must be an 
even number of bars left, half of type n and half of type U. As the bars n have factor 
g-OT(±fl) Ij factor we have full cancellation again and have total 

factor 1. 

For the remaining case cr e 2^*^j,„((x)) \ S'^'^(w) and (-l)lh->’ll = _i and to e E[xo —yo], 
we have xq and yo in different sublattices. We can see as last time that the factors from 
the ‘extra bars’ (that arise from each half of the loop having a different number of bars) 
will cancel as again there are equal numbers of n and U. For the remaining bars there 
are an odd number in each half of the loop, this means we can make a matching for all 
but two of the bars. The factors from bars in the matching will cancel each other. For 
the remaining two bars one is a n with factor and one is a U with factor 
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(the sign of a is opposite due to the spin flip at xo at time 0). This means the overall 
factor is (±0^ = -1- This completes the proof. □ 


In light of Proposition |5. 1 1 the following proposition can be proved in the same way as 
Theorem 3.2 in Il23l . 


Proposition 5.2. The partition functions Z' 


( 1 ) 

A./? 


are given by 




(5.4) 


We also have the following identity 


Tr(cr^ ® l);,(cr" ® l)ye'^"A 


-m' - 




P(dw) ^ 


Y\ Zxi,yj(o-,OJ) 


\{xi,yj}e6 


^ An ^ V< 


(5.5) 


where cr^. is the value of a space time configuration, cr, at time 0 and site Zi- 

With the important details understood we can calculate some correlations in terms of 
probabilities in the loop model. The most important correlations here are the Neel and 
nematic correlations (Proposition |5. 3 1 (a) and (b) respectively). 


Proposition 5.3. For i, j — x + y, i + j and A bipartite 

(a) = (-l)ll-^IIP(x(^yQ), 

(b) + 5P(THHj + 5 p(!Z!°) + 

(c) {S[sis‘ySi)Afi = i [-(-l)ll-'-yllP(x,^yo) + P('“£“), 

(d) = I [(-l)ll-'-''IIP(x(^yQ) +P(]“£“), 

(e) ^+i [p(]v)+p(? j;) - p(;;2;;) 


Proof We will calculate the correlations in order. First note that each S‘ plays an 
equivalent role, hence cyclic permutations of the indices (1,2,3) does not alter the 
expectation. Using this together with {S‘S^Y - ±{S^S‘) (the sign depending on the 
value of i and j) means we can take i = 3 and j - I- For each we will expand using 
(|A8) and ( |4T2l l. 

Proof of (a). First 

{S lSpA,g = ®l®cr^®l+cr^®l®l®cr^ + l®cr^®cr^®l + l®cr^®l ® 

(5.6) 

We see that due to sites zo and Zi being interchangeable for z e A each of the four terms 
in the sum have the same expectation. We also know from Proposition |5^ 

Tr{o-^ l)Yo-^ ® l)ye~^^^ ^ f p(daj) ^ cr^oO-y„. (5.7) 
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We note that the integral differs from zero only on the set where xq and yo are connected. 
If X and y are in different sublattices the product of spin configuration values is if 
in the same sublattice the product is We can deduce that 

(SlSpAfi = (-l)""-^"P(x,^yQ). (5.8) 


Proof of (b). For the second correlation 


(Rif = (cr^ ® 1 + 1 ® 0-3)2 = I -1 ® 1 + 2cr" ® cr^ 


(5.9) 


We see that expanding as before gives 




f Z l^+2cr. 


xo^xi I “ P(*^0—-^l)- 


(5.10) 


From this and the fact that {{Slf)^^ - 5 (S;c ■ Sj:)a^ = | we can deduce that 


P(jco—Jci) = 


(5.11) 


For the hrst term in the correlation we again note that {(S lf{S lf)Afi - {(Riyi^f/Ajs- 
We then calculate as before: 


(RiyiRiy =(0-3 ® 1 + 1 ® Cr3)2(o-3 ® 1 + 1 ® cPf 


1 


= -1 ® 1 + 2o- ® cr -1 ® 1 + 2o-^ ® cr 


1 


(5.12) 


- 1^-1®'3 +0-3®0-3®l®l + l®l®0-3®0-3+ 4(o-3)®^ 
Now following through the same expansion as before we have 


<( 52 ) 2 ( 53)2 


)a,/j J Z 




^ -I- CTxoCTxj + CTyo'^ri ^O-xoO-xjO-y^O-j,! 


(5.13) 

Using ( |5.11| l and noting that the last term in the sum requires either two loops con¬ 
taining two of the sites xq, x\,yQ,yi each or one loop containing all four sites to give a 
non-zero contribution to the sum overall (if one site is not connected to any other its 
spin value can be + j independently of other sites, averaging the integral on this set to 
zero) we have 


<(53)2(53)2)^^ - {(Slf)Afi{{S]f)Afi = + ^P(jV) + + ^KM)- 

_ . -M . >. 

The probability P^ " ={= comes with twice the weight because there are two ways to 

connect both sites at x to different sites at y (but only one way both sites at x can be 
connected and both sites at y can be connected). 

Proof of (c). For the third correlation we use the same expansion 
(SqS gS1 }a./3 - -^7’r(cr'cr3®lH-cr*®cr3)j^(cr'o-3®lH-cr'®cr3)j,FAe~^^A ■ (5.15) 
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The factor 4 has come from grouping together terms such as cr* ® cr^ and cr^ ® cr* that 
have the same expectation. A useful observation at this stage is that cr'cr^ = 
Calculating further and noting that the two cross terms in the above product have the 
same expectation we see 


/ 1 

{S qS Is '5 = 4(--cr^®l®cr^®l - icr^ ® 1 ® cr' ® cr^ + cr' ® cr^ ® cr^ ® cr^ \ 

\ 4 lAfi 

(5.16) 

From the symmetric roles of cr' for i = 1,2,3 and part a) we know the first term is 
_ (-iy' ^" p(.xo—yo)- For the second term we need {cr^ ® 1 ® cr' ® cr^)^*^- This is the 

expectation of a matrix with purely imaginary entries, due to the one appearance of cr^. 
Now we note three pieces of information that allow us to calculate this expectation. 
All the matrices ,P/^,cr^,cr^ and cr^ are Hermitian. The matrices cr' are acting 
on different sites in A and hence they commute, e a and Pa commute and have real 
entries. This means taking the adjoint of the operator leaves the expectation unchanged. 
Because the operator is purely imaginary we should obtain the negative of what we 
started with on taking the adjoint. Hence the correlation must be zero. 

For the last term we expand as in Proposition |5 . 2| and obtain 


(cr' ®cr^®cr' 


-(Ij- I P(dw) ^ Z;c„,j,„(cr,w)<cr._o+|cr'®cr3®o-'®0-V-,0-> 

Jot'll 

(5.17) 


Here cr. o± denotes the full spin configuration for some cr € 'Lx„^y„(io) at time 0+ re¬ 
spectively. Also note that as cr* flips spins and cr^ does not the set of space-time spin 
configurations 2^2>^^(cu) is the correct set. We could expand the set of configurations 
we sum over to include configurations that flip spin at sites xi and y\ at time zero but 
these would not be compatible with cr^ acting at time zero at those sites hence they 
would not contribute. Recall that a loop that contains a site that spin flips at time zero 
cannot contain only one such site, hence the set of configurations that contribute to the 
integral is E{xq —yol- Again the set of configurations where one of the sites xj or yi is 
not connected to any of the other three does not contribute to the integral. Combining 
these two facts we see that the only sets of configurations that contribute to the integral 
are those where there are two loops each containing two sites (one with xq and yo and 
the other with xi and yi), or one loop containing all four sites. For the case of two 
loops there is one factor of Zxa,yo(cr, oj) - (-l)ll'"''ll from the loop containing xq and yo 
(where cr* acts). Another factor of (-l)lh-3'll comes from the loop containing xi and yi 
and the condition that the spin flips on crossing a bar. Note that for the first loop there 
is no such factor coming from spin flips at bars because cr*| + i) = + 5 I + \) hence 
there is a factor of H- 5 regardless of the spin value at the site. For the case of one loop 
containing all sites the order that sites occur in the loop is important; this is because 
both cr* and cr^ are acting at sites in the loop. If, when following the loop, the site yi 
appears directly before or after the site x\ then the section of loop between these sites 
follows the normal rule of flipping spins at bars (or if we follow the loop the other way 
we pass through two spin flips at time zero as well, these cancel each other out as far 
as the product of spins at sites xi and yi is concerned). This means we have a factor 
of (-l)lh-rll as before. If one of the sites xq or yo appears between sites xi and yi on 
the loop the effect of the extra spin flip changes the sign of the factor coming from 
the product of spins, giving a factor of -(-1)II’^~^IL As before we also have the factor 
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(cr,w) = (-1)"^ in both cases. This means the correlation is 


^xa,ya 


{cr^ ® cr^ ® cr' ® “ 



(5.18) 


Recall that the arrows in the events show the direction that the loop is traversed. From 
this we can finally deduce that 




-l)ll"-yllp(x,^yQ) + + p(?^) - p(|x|) 


(5.19) 


Now we note that the last two probabilities are equal (swap yg yi). 

The correlations (d) and (e) follow easily using the same techniques and considerations 
as above. □ 


From this we can easily obtain some bounds on these correlations that are potentially 
very difficult to obtain without the loop model. 


Corollary 5.4. For i, j — 1,2,3, x + y, i + j and A bipartite 

(a) <(5')2(5')2 )a,^ - <Ti + |(-l)""-""<5i5')A^ 

(b) {S‘,sis‘ySl)A,/t < 5((-l)""-"" - l)<5iA')A,/, 

(c) {S‘^SiS^ySpAj3 < i((-l)"-'-"" + IX^i^pA./^ 


(d) {S‘,sis‘ySi)Ajtl^ 

(e) <5i.5i5^5pA,/?| 


~ yll is odd 
“ yll is even 

if 11-^ “ yll is odd 
if 11^ “ yll is even 


Proof. All inequalities are immediate from Proposition 


5.3 when we note that E 

1 ° T° 

H^l 

LJ 

XI y, 


is a sub-event of ^[xo—xj and E 


is a sub-event of ^[xo—yd. 


□ 


Other inequalities of interest involve correlations between nearest neighbour points. 
Equation (29) in ll20l allows us to obtain the following bound in the ground state {f 

oo) 

1 272 - 37i 

Now looking at Proposition |5.3| (b) for ||x - y|| = 1 (say x = 0, y = ei) we see that if 
7i = 0 then the event P(Oo— ei„) puts us into the case of one of the last two probabilities. 
Ignoring the first probability (as it is difficult to control) we obtain (for 72 > 0 = 7i) 

- {{Slf)Ad(Slf)A,(t ^ (5-21) 

This bound is positive for 7 < 4, however it was not sufficient to deduce nematic order 
from a theorem analogous to|6.1|but concerning the nematic correlation function. 
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6 Occurrence of Neel Order 


6.1 Setting and results 


We take the cubic lattice in with side length L, denoted A/,, with periodic boundary 
conditions. The edge set, will consist of pairs of nearest neighbour lattice points. 
Precisely 


Al = 




c Ail \\x - yll = 1 or jx,- - y,| = L - 1 for some / = 1,..., d]. 


( 6 . 1 ) 

( 6 . 2 ) 


For the main theorem we need to introduce two integrals, they come about due to 
similar considerations as in ina 


Kd^ 


— f -f 

f 


(In)' 

1 

(^ 


d \ 

COS ki 


j s(k + 7 t) 
' E(k) 


Ak, 


\e{k + ji) 
' e(k) 


dk. 


(6.3) 

(6.4) 


Here (■)+ denotes the positive part and £(k) - 2 - cos(ki)). 


Theorem 6.1. Let d > 3 and J\ <Q < J 2 , for L even we have the two bounds 


lim lim - 

/3~»oo L-.00 |Ai| 


xehi 


> < 


VP(Oo—eio)' (VP(Oo—eio)' - h 



(6.5) 


Positivity of this lower bound implies Neel order for those values of J\ and J 2 in the 
spin-1 system at low enough temperature. This result is also equivalent to the occur¬ 
rence of macroscopic loops in the loop model. 

Of course we see that for -J\,J 2 > 0 the positivity of the lower bound doesn’t depend 
on the value of + J^, only on the ratio -7i/72- This means there corresponds an 
angle, measured from the Ji axis, such that for angles less than this we have proved the 
existence of Neel order/macroscopic loops. The bound is positive if 


-ei„) 


Kd^ -J1+4J2 


Vp(Oo— ei„) > Id 



■h 

Ji ■ 


( 6 . 6 ) 


One of these is certainly satisfied if IdKd < (-47i)/(-7i + 4J2). A table of values of 
Id and Kd for various d is presented in ||2^ . \f J^ + - 1 this is the case in = 3 for 

Ji < -0.42, d - 4 for J\ < -0.28 and = 5 for 7] < -0.22. 


A similar theorem concerning nematic order (corresponding to correlation (b) in 5.3 1 
can be proved using the same methods. Unfortunately showing that one of the lower 
bounds obtained was positive proved difficult due to the seemingly unavoidable issue 
of bounding more complicated connection probabilities from below. 
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6.2 Proof of Theorem 16.11 


The result can be proved directly using the loop model. However in light of the proof of 
reflection positivity for the J 2 interaction m, we can complete the proof more directly, 
appealing to Proposition |5 .3 1 when required. To begin we define a Hamiltonian 


= Z hjl. (6.7) 

\x,y}e£, xeAi 


denote its Gibbs states {■)'. Note that for h = 0 this Hamiltonian is unitarily equivalent 
to as can be seen from defining U - nxeAj, where c A/, is the odd 

sublattice. It has been shown Eng that this Hamiltonian is reflection positive for 
7i < 0 < J 2 - For external field v = (vt)xeA£ for V;c 6 R define 


= ^a;,o + -^1 Z 

xeAi 


and 

Z(v) = (6.9) 

By the usual methods for reflection positivity, using m to handle the J 2 terms we 
can show that for any v we have Z(v) < Z(0)e~'^‘^^''’^''\ By choosing v according to 
Vx - rjcosik ■ x) for rj small we can recover the infrared bound in ||8l : 




1 


i-2Ji)sik) 


( 6 . 10 ) 


where 

{A,B)Duh^^^\ AsTrA*e~^^^L , ( 6 . 11 ) 

is the Duhamel inner product. Now we use the Falk-Bruch inequality 

1(A*A + < 1 ^(A,A)ouh + ^(A,A)Duh. (6.12) 


After some calculation and use of the SU(2) symmetry we find 


''ZV ■ 0/ t 


= \K:^\e{k + n)[2J,{Slsiy^';’^ 

-mslfiS^yyy - 8(5^535^,53^)1';^^ -4(515353^5],)^;^^)). 

(6.13) 

Using Proposition |5.3| we can obtain the bound 

{[S\,[Hy\Sl]\yy < |Ade(k + 7r)(-27i +872)P(0,^xo ). 

Now we use the Fourier identity 

1 


(6.14) 


IAlI 


^(- 1 )«( 535 ])~ = ( 535 ])~-^ ^ e‘’^y{Syiy{k) ( 6 . 15 ) 

xeA, ' teAniOl 


' teAIMO) 

with y = 0 and y - ei. We see from using unitary operator U with ( |6.10| l and ( |6.14| i 
have after taking limits the bounds in Theorem 6 .1 


we 
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